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Abstract 

We study the conditions for neutrino oscillations in a field theoretical approach by 
taking into account that only the neutrino production and detection processes, which are 
localized in space around the coordinates xp and xd, respectively, can be manipulated. 
In this sense the neutrinos whose oscillations are investigated appear as virtual lines con- 
necting production with detection in the total Feynman graph and all neutrino fields or 
states to be found in the discussion are mass eigenfields or eigenstates. We perform a 
thorough examination of the integral over the spatial components of the inner neutrino 
momentum and show that in the asymptotic limit L = \x*d — xp\ oo the virtual neutri- 
nos become "real" and under certain conditions the usual picture of neutrino oscillations 
emerges without ambiguities. 



1 Introduction 



In particle physics, oscillation phenomena have so far only been observed in the K°K° 
and B°B° systems. Analogous oscillations are expected to take place between neutrinos 
of different lepton flavours |l| if neutrinos are massive particles and if individual lepton 
numbers are not conserved. This non-conservation manifests itself in a neutrino mass 
matrix which is non-diagonal in the weak basis, i.e. the basis where the mass matrix 
of the charged leptons is diagonal and the charged current interactions of the standard 
model (SM) is of the form 

£cc = -^^j: u iL r£ L w++h.c^j. (i) 

Though neutrino masses and mixings rank among physics beyond the SM it is nevertheless 
reasonable to assume that neutrino production and detection are within the province of 
the SM interactions. In this sense the left-handed neutrino fields in eq. (HD are linear 
combinations of mass eigenfields, i.e. 

u eL (x) = U ej v jL (x) (2) 

with a unitary 3x3 matrix U and neutrino fields Vj with mass rrij. 

The standard approach to neutrino oscillations (see, e.g., ref. 0) is important for 
its simplicity and its physical insight but is not sufficient for a complete understanding. 
(See, e.g., ref. || for a list of its conceptual problems.) It is well known that it works only 
in the extreme relativistic limit. Wave-packet treatments |4]] are not totally satisfactory 
either without knowing size and form of the wave packets. However, it is clear that, if 
the particles involved in the neutrino production process are assumed to have definite 
4-momenta, the neutrino is forced to have a definite 4-momentum too and hence is 
in one of the mass eigenstates. Therefore, in order to observe neutrino oscillations one 
needs a sufficient spread of momentum (or energy) of the particles in the production and 
detection mechanism. It has been stressed in the literature |3|, |5| that the correct way to 
avoid all difficulties and ambiguities associated with neutrino oscillations is to concentrate 
on those things which can really be manipulated or observed like the particles responsible 
for neutrino production and the target responsible for neutrino detection. This amounts 
to including the production and detection mechanisms localized at the coordinates xp 
and xc, respectively, in the calculation such that the neutrinos propagating from xp 
to xp, are virtual particles or intermediate states (see ref. || for a field-theoretical 
treatment and ref. 0] for a quantum-mechanical treatment along these lines). 

In this paper we choose the field-theoretical approach and apply it to the antineutrino 
production as it occurs in a reactor and assume that the antineutrinos are detected by 
elastic electron scattering. (This is, e.g., the situation in the MUNU experiment 0.) The 
antineutrino flux is produced by the /3-decay of the fission products in the reactor. We 
simplify this mechanism by taking only the /3-decay of the neutron and assuming that 
both neutron and proton are bound to nuclei at rest. This merely serves to keep our 
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notation as simple as possible and has no impact on the essence of our discussion. The 
main points in our treatment of neutrino oscillations are the following: 

• The neutron and the proton are bound to nuclei in stationary states localized at 
the coordinates xp whereas the target electron is bound in an atom at the point 
xp>. Neutrino production and detection are macroscopically separated by a distance 
L = \L\ with L = — xp. 

• All particles in the final state are described by plane waves. This is justified for the 
electron in the detector and all the more for the other electron from the /3-decay 
and the antineutrino which are both not observed. In the event rate it has to be 
summed (integrated) over a complete set of states for these unobserved particles 
and plane waves obviously serve this purpose. (This is one of the points which 
distinguishes our work from ref. H.) 

• The internal (anti) neutrino is represented by its Feynman propagator. The integra- 
tion over this internal line with momentum q is coupled to the macroscopic distance 
between neutrino production and detection by / d 3 q exp(—iq- L) (q 2 — m 2 + ie)^ 1 . . . 
for a neutrino with mass rrij. We use a mathematical theorem which is proved in 
the appendix to obtain the leading term of the amplitude in the limit L — * oo. We 
will see that in this asymptotic limit the antineutrino is "forced" to be on its mass 
shell and the momentum q is aligned with —L. 



2 The oscillation amplitude 

We consider the following weak process 

n — > p + e~ + v 

\ 

^ + e D ~ ¥ & + e D (3) 

occurring through the intermediate propagation of an antineutrino, where n (neutron) and 
ejj (electron) are the initial particles. As already mentioned we assume that the neutron, 
the proton and the target electron are in stationary states. Since neutrino production 
and detection are localized at xp and xp,, respectively, the spinors of the initial particles 
and the proton can be written in position space as 

ip p (x) = ip p (x — xp) exp(—iE p t), 
if) n (x) = ij) n {x- xp) exp(-iE n t), 
^eD{x) = ipeD{x - x D )exp(-iE eD t), (4) 

respectively. The functions ip p , n ,eD(y) are peaked at y — 0. The final particles will be 
described by plane waves. 
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The weak interaction Lagrangians relevant for production and detection are 

Q 

C P (x) = -^cos$ c J2 U ej e(x)7 A (l - 7 5 )z/ j (x)p(x)7 A (l - g^)n(x) + h.c, 

v 2 j 

£ D (x) = E {e(x)^P L e(x) (2U* ej U ek + 5 jk (g v + g A )) + 

V 1 j,k 

+e(x)^P R e(x) 5 jk {g v - g A )j z^(x)7 M (l - 7 5 )z/ fc (x), (5) 

respectively, where Gf is the Fermi coupling constant, $ c the Cabibbo angle, gy = 
2sm 2 9 w - 1/2, g A = -1/2, W the Weinberg angle, g « 1.26 and P fijL = (1 ±7 5 )/2. All 
fields occurring in eq. (|5|) are mass eigenfields. The amplitude for the process (^|) with an 
antineutrino of mass m k in the final state is given by 

A k = (p,u k (p' u ),e~(p' e ),e]y(peD)\T[J d 4 xi J d 4 x 2 C P (x 1 )C D (x 2 )]\n, e D ) . (6) 

With the neutrino propagators of the mass eigenstate neutrinos 

we obtain the amplitude 

A k = ^|^£, E J J d * X2 J 0_ e -M^ 2 ) 

exp {-i{E n - E p )t x - iE eD t 2 ) exp (ip' e ■ x x + i(p' u + p' eD ) ■ x 2 ) 

j x {x x - x P ) uMh x a - i 5 )u ej ( V(i - tbK(^) 

\q z — rrij + «e/ 

UeifeDhfj, Pl (2U* ej U ek + 5,-fc(#y + 9a)) + PR$jk(gv ~ 9a) ^eD{x 2 - X D ) (8) 

where 

J x {xi - fp) = ^ p (fi - £p)7a(1 - gib)^n{xi - x P ) . (9) 

Note that the amplitude (§) is not antisymmetric with respect to the exchange of the 
final electrons because the electron generated in the /?-decay never reaches the detector 
or, in other words, it is absorbed in the reactor and thus takes part in the formation of 
final states orthogonal to single electron states. 

The integration over t\ and t 2 in eq. (||) can easily be carried out leading to the delta 
functions 8{q Q + E 2 ) and S(q + Ei). Furthermore, we use the relation 

dxe~ ikx f(x + b) =e lkb f(k) , (10) 
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where / is the Fourier transform of f, in the integrations over x\ and x 2 . Thus we obtain 
Ak = i^H^ £ / A S(q + E 1 )5(q + E 2 ) 



2(2tt) s 

Upi ~ 0) ^e(pi)7 A (l - ls)U ej ( / + ™f • W - 76KW 



u e (p'eD)Y \P L (2U* ej U ek + 5 jk (g v + g A )) + P R S ]k (g v - 9a)\ $ eD (p2 + g) (11) 



(12) 



where we have defined 

E x = E n -E p - E' e , p x = p' e , 
E 2 = E' U + E' eD - E eD , p 2 = p' eD + f u . 

Now only the integration over q remains. Since our initial and final states are all energy 
eigenstates, the integration over q° leads to the usual 5-function 5(E 1 — E 2 ) expressing 
that the initial energy is equal to the final energy of the total process (|3|). The analogous 
5-functions corresponding to momentum are smeared out by the initial momentum dis- 
tributions and therefore we have a non-trivial d 3 q integration in the amplitude. We can, 
however, take advantage of the fact that this integration actually amounts to calculat- 
ing a Fourier transform and evaluating it at the coordinate L, the macroscopic distance 
between neutrino production and detection (see eqs. (f|)). Hence we can apply the theo- 
rem proved in the appendix (see eqs. (A.l) and (A. 2)) which enables us to calculate the 
leading term of the amplitude for large L. 

Let us for a moment dwell upon the mathematical requirements of this theorem in the 
context of the phyiscal problem under discussion. The wave functions eq. (|j) represent 
bound states. It is therefore reasonable to assume that they fall off exponentially and 
consequently the Fourier transforms J\ and ip e D are infinitely many times differentiable. 
This meets the requirement of differentiability in the theorem. It remains to discuss the 
behaviour of J\ and ip e n for \q\ — > oo. If bound states in coordinate space are also 
infinitely many times differentiable in addition to their exponential fall-off then their 
Fourier transforms have the same properties. However, bound state wave functions are 
in general not everywhere infinitely many times differentiable in coordinate space. A 
typical example is given by the ground state wave function of the hydrogen atom which 
is not differentiable at the origin. Its Fourier transform nevertheless decreases with the 
fourth power of \q\ which meets the conditions for the decrease of the function in the 
theorem. Therefore the physics under consideration seems to comply naturally with the 
mathematical requirements. 

Noting that the constant A of the theorem corresponds to q$ — rrij, eq. (A. 2) tells us 
that the asymptotic limit L — » oo of eq. ( JTT| ) is given by 

(l 2 rrmi9 1 

A? = -i c ~ 5(E 1 - E 2 ) e -m-z P -iP2-x» 

Y,e iqjL J\(Pi + Qjl) u e (Pi)l X (l ~ 75)£4;(*)7° + ?/'7 + ^i)7 M ( 1 ~ IsWipl) 
j 

u e {p'eD)Y \Pl (2U* ej U ek + 5 jk (g v + g A )j + P R 5 jk (g v - g A )} i) eD (p 2 - q 5 l) (13) 
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with 

/ = -i-, q = -E 1 = -E 2 and <?j = y/qfi ~ mj . 

This equation can be interpreted in such a way that .41° is a sum over terms which 
contain a real antineutrino of mass rrij. The reason is that the 4-vector pj defined by 

p°.=E u = -q = E 2 = E' eD + E' u -E eD >0 and = g/ (14) 

can be identified as the 4-momentum of this antineutrino. The spatial part of pj correctly 
points from xp to Xr>. Furthermore, with a suitable renormalization of the 4-spinors v 
the identity 

go7° + Qjl ■ 7 + m j = ~ i>3 + m i = - £ > s )%(Pj> s ) ( 15 ) 

shows that apart from the factors 1/L and exp(iqjL) the amplitude ^4^° is just the sum 
over the products of production and detection amplitude of antineutrinos with mass rrij. 

3 Oscillating amplitude and cross section 

We expect that under certain conditions the cross section for the reaction @ (n + — > 
er + + P) exhibits an oscillatory dependence on the distance L. Considering first the 
asymptotic amplitude A™ and assuming that E v ^> rn^Vj we can make the following 
remarks and observations to this effect. (Most of them have already been discussed in 
the context of the quantum-mechanical approach of ref. 0.) 

• In order to obtain ^4°° we have performed the asymptotic limit L — > oo. Looking 
at theorem proved in the appendix it is evident that "L large" means E u L/%c ~ 
E v L/2 ■ 1CT 13 MeV • m ^> 1 where E v is an average antineutrino energy and L 
is measured in meters. For every thinkable neutrino experiment this is very well 
fulfilled and corrections to Af are suppressed by (E u L/hc)^ 1 ^ 2 . 

• The factor 1/L in the asymptotic amplitude corresponds to the geometrical decrease 
of the neutrino flux by 1/L 2 in the cross section. 

• Looking at eq. ( jl3D we conclude that neutrino oscillations with masses m^m^ can 
only take place if 

I Qj ~ Qk | ^cr,/ and | qj - q k \ < a D 

where oj and op are the widths of the functions J and if) e r>, respectively. In coor- 
dinate space this simply means that the widths of J and ip e £> must both be smaller 
than the oscillation length Z osc = A-KE v /\m 2 j - m|| « 2.48 m • ^(MeV)/Am 2 (eV 2 ). 

• To evaluate this condition for reactor neutrinos we note that hc/aj> 10~ 15 m and 
/ic/<7d^10 -10 hi. From the latter estimate one gets the condition Am 2 <(100 keV) 2 

I- 
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• The coherence length is infinite in our calculation since all initial and final states 
are stationary. 

• Our calculation was based on the assumption that neutrinos are of Dirac nature. 
In the case of Majorana neutrinos the diagonal vector currents in £,£> eq. (|5j) are 
identically zero. As a consequence, in Ak eq. (H) the second 1 —75 has to be replaced 
by —275 for j = k. This introduces the well known difference of order mj/E v 
between the amplitudes of Majorana and Dirac neutrinos ||. 

Finally we come to the cross section or event rate in order to make contact with 
the usual oscillation picture. If \qj — qk\ <C min(crj, ctd) is fulfilled V j, k in addition 
to the assumption that all neutrinos are ultrarelativistic, then we can take the limit 
rrij — > Vj in all terms of A^ except the exponential factors exp(iqjL). Then calculating 
the cross section one would arrive at the same result as can be obtained by the following 
heuristic consideration: The probability to find a neutrino at a distance L from the 
source is given by Pi = | J2j U e jUlje tqjL \ 2 . Thus the number of events in elastic ve~ 
scattering at a certain neutrino energy E u is proportional to J2e Pt <J {^t e ~\ E v ) where the 
a(i>ie~; E u ) (£ = e, /i, r) are the elastic scattering cross sections as given by the SM. Note 
that E,. can be reconstructed via 




E v = m e 4/1 H — — cos a — 1 (16) 



where T = E' eD — m e and a is the angle between L and PeD, the momentum of the recoil 
electron. Thus the event rate is given by 

"' A ^ {| £ l^f^r (^ e ~; E „) - afre-; E v j) + a(^e-; E v )} (17) 



dE v 



,1 



with a normalization constant Nq. 

In conclusion we want to stress once more the central importance of the theorem 
demonstrated in the appendix. It is not only relevant in the context discussed here but it 
shows in general under which circumstances the total amplitude for particle production 
and its subsequent scattering factorizes into a product (or sum of products) of production 
amplitude and scattering amplitude. In the case of neutrino oscillations it leads to a 
transparent and simple field-theoretical treatment without ambiguities or conceptional 
problems 0: all neutrino states or fields are mass eigenstates or eigenfields and the 
neutrino "wave packet" is totally determined by the production mechanism. The method 
described here also provides the correct answer in cases where the standard approach to 
neutrino oscillations fails. We want to mention, however, that some more work has to be 
done to extend our method because it is not straightforwardly transferable to the case of 
accelerator neutrinos and it has the deficiency of an unrealistic infinite coherence length 
for neutrino oscillations. 
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Appendix 



Theorem: Let $ : R 3 — > R 3 be a three times continuously differentiable function such 
that $ itself and all its first and second derivatives decrease at least like 1/q 2 for \q\ — > oo, 
A a real number and 

J(L) = J d 3 qme- i ^ A _l 2+ . e (A.1) 
Then in the asymptotic limit L — \L\ — > oo one obtains for A > 

J{L) = e iyaL + C(L- 3/2 ) (A.2) 

whereas for A < the integral decreases like L~ 2 . 

Remark: In order to make the proof of this theorem transparent we will first introduce 
three lemmas and then divide the proof into several steps. 



Lemma 1: Let / be a three times continuously differentiable function and 

J(r) = I* d6 sin 6 e~ ircos 9 f{0). (A.3) 
Jo 

Then in the asymptotic limit r — > oo one obtains 



J(r) = _*-' r /(0)-e'7M +orn (A4) 
ir 



Proof: We first perform a partial integration resulting in 



/(r) = _e-'7(0)-e"7(.) _ i/i(r) (A g) 



with 

- ir cos 9 



U r ) = / d9e- trcos0 f(6). (A.6) 

JO 

It is convenient to split hij) into the two parts 

Ur)= [ W/2 d6e- ircosd f (6) + ["^ d6 e ircos 9 f'(ir - 6). (A.7) 
J o io 

In the following we only discuss the first integral in this equation since the second one is 
treated analogously. Defining g{6) = (f'(0) — f'{0))/ sin 6* we obtain 

d9e~ irco8e f(e) = [ W/2 d6e- ircos9 f'(0) + ' d6 sinO e- ircos9 g(6). (A.8) 
Jo Jo Jo 

It is easy to show that g is twice continuously differentiable in the intervall [0, 7r/2] as a 
consequence of the three times continuous differentiability of /. Therefore one can apply 
a partial integration to the second integral on the right-hand side of eq. (|A.8|) in the same 
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way as to eq. 



and thus this integral decreases like 1/r in the asymptotic limit (see 
eq. A.5|) . Together with the 1/r in front of h{r) in eq. (|A.5|) there is an overall decrease 
of 1/r 2 . Therefore it remains to consider the asymptotic behaviour of 



tt/2 



d6e 



dz 



l+w/r 



dz 



2 \J0 VI - ^ Jn/r _ {z _ 2)2 

The right-hand side of this equation allows to deduce the upper bound 



7I-/2 



d6e 



< 



1 

< - 
~ 2 



Tr/r 



dz 



VT 



\ i-l+w/r 

+ / dz 



n\2 



1 

- / dz 

Jn/r 



arcsm 



1 - u - - 



n\2 



'27T 7T 2 



2 ' 



(A.9) 



(A.10) 



Thus the integral of eq. ( |A.9| ) is bounded by a function with asymptotic behaviour 
0(1/ y/r) and lemma 1 follows. □ 



Lemma 2: Let / be a three times continuously different iable function and w be a real 
number with \w\ < 1. Then the integral 



IJr)= / d9 sm9e- ircos9 sm(wrcos6)f(6) 
Jo 



has the asymptotic behaviour 

w cos wr + i sin wr e ir 



Ur) 



-/(0)- 



1 - 



, . w cos wr — i sin wr 

— / W ; 2^ + 0(r- 

r 1 — wr 



3/2n 



(All) 



(A.12) 



in the limit r — > oo. 



Proof: The proof of this lemma follows from lemma 1 by using sin a = (e m — e ta )/2i 
for a = wr cos 8. □ 



Lemma 3: Let ip : R 3 — > R 3 be a twice continuously differentiable function such that 
the function itself and its first and second derivatives are absolutely integrable. Then the 
integral 

I{q) = J d 3 q4j{q)e-^ Z (A.13) 

decreases like L~ 2 for L — > oo. 
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Proof: According to the assumption, is absolutely integrable and with a twofold 
partial integration we obtain 

J d 3 q (Aif>)(q) e-^ £ = —L 2 J d 3 q^{q) e~^ £ (AAA) 

which proves lemma 3. □ 

Proof of the theorem: The main point in the proof is to meticulously take care to 
perform the limit e — > before investigating the asymptotic limit L — > oo. 

Step 1: For A < we apply lemma 3 and thus obtain the second part of the theorem. 
In the following A will always be positive. 

Step 2: Next we split J(L) into the integrals 

Ji(L) = -jVAj s2 dn<S>(VAn)e- iVJnZ (A.15) 
originating in the S— function of the limit e — > and 

ML) = ~J ^me-* \ q2 q *~f +€2 (A.16) 

which represents the principal value of this limit. S 2 denotes the 2-dimensional unit 
sphere, q = \q\ and n = q/q. 



Step 3: We choose two numbers 5 and r) such that 0<5<rj<\/A and a symmetric 
C°° function h such that < h(v) < 1 for all v G R, h(v) = 1 for \v\ < 5 and h(v) = 
for \v I > 77. Then J 2 = Sfe=i <^2fc with 

J 21 (L) -/ d 3 qe-^(m-H^An)h(q-VA)) ^/'y^^ (A.17) 

J 22 (L) = - J^dtt J™ dqe- l * Z <5>{VAn)h(q-\/A) ■ 
( q\q 2 ~A) 2A 3 ' 2 {q-^A) \ 

J 23 (L) ^ - / dfl T dge-^ *(VAn) h(q - VX) if^Z^ 2 ■ (A.19) 
JS 2 7o 4A(q - VA) 2 + e 2 

In these three integrals the limit e — > can be performed. In J 21 one obtains an integrand 
which is twice continuously differentiable and therefore J 2i decreases like L~ 2 for L — > 00 
according to lemma 3. The function within the parentheses of J 22 gives (q + y/A/2) /(q + 
\/A) in the limit e — > and with the properties of h the same asymptotic behaviour as 
that of J 2 i results. Finally, with the transformation v = q — \fA the integral J 23 can be 
cast into the form 

J 23 (L) = 1 -^aI dn $(VAn) e- iVJnrL f dv h(v) ' L) . (A.20) 

2 J s 2 J—r) v 
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Step 4: To perform the integration over S* 2 we take 



L = L 



( \ / sin 9 cos 

= Le z and n(6,(f))= sin 9 sin (p 
^ 1 / \ cos# 



Then the application of lemma 1 with r = \H\L shows that 



(A.21) 



J\{L) = {e~ iVjL ^{VAe z ) - e iVjL <f>(-VAe z )) + 0(L- 3/2 ). (A.22) 



Step 5: By virtue of lemma 2 and taking w = v/y/A, the asymptotic behaviour of J 2 3 
is given by 



ml) = % r dv- ^ 



e - iVjL ^(^Ae z ) l-^=cos(vL) +ism(vL)j - 



-e iVjL $(-VAe z ) l-^=cos{vL) -ism(vL) 



+ C(L- 3 / 2 ). (A.23) 



The terms with cos(vL) decrease faster than any power of 1/L by a corollary of lemma 
3 whereas the 5— function property of sm(vL)/nv leads to 



7T" 



J 23 (L) = -— (e^ L $(VAe z ) + e iVjL $(-VAe z )) + 0(L~^ 2 ). 
The correction to the 5— function is proportional to 



(A.24) 



sin(uL) 



(A.25) 



and therefore decreases faster than any power of 1/L as can easily be seen by repeated 
partial integration. Finally, the theorem is obtained by summing the results of steps 4 
and 5. □ 
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